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A theory of collective states in a magnetically quantized two-dimensional electron 
gas (2DEG) with half-filled Landau level (quantized Hall ferromagnet) in the pres- 

m ■ 

ence of magnetic 3d impurities is developed. The spectrum of bound and delocalized 



spin-excitons as well as the renormalization of Zeeman splitting of the impurity 3d 
levels due to the indirect exchange interaction with the 2DEG are studied for the 



specific case of n-type GaAs doped with Mn where the Lande g-factors of impurity 



and 2DEG have opposite signs. If the sign of the 2DEG g-factor is changed due to 
external influences, then impurity related transitions to new ground state phases, 
presenting various spin-flip and skyrmion-like textures, are possible. Conditions for 
existence of these phases are discussed. 
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I. INTRODUCTION 

In a strong magnetic field the two-dimensional electron states in semiconductor 
heterostructures 1 transform into Landau states with a completely discrete energy spec- 
trum. This diamagnetically quantized two-dimensional electron gas (2DEG) possesses many 
remarkable features including Quantum Hall effect. 2 The role of impurities in the thermo- 
dynamic, optical and transport properties of 2DEG is extremely important. Among many 
facets of this problem we choose for discussion in this paper the formation of impurity related 
collective excitations in a magnetically doped quantized 2DEG in the case of odd integer 
filling factor v — 2n+l. In a pristine state 2DEG with odd v is in a Quantuzed Hall Fer- 
romagnet (QHF) regime with nondegenerate ground state characterized by the total spin 
quantum number S = and maximum spin projection S z = S. (N^ is the magnetic- 

flux-quanta number.) Different branches of the excitons are well distinguishable among the 
low-energy excitations. They are classified as spin waves (spin excitons), magnetoplasmons 
or multi-exciton states depending on the spin and orbital quantum numbers. 3-10 Besides, 
low-lying collective half-integer-spin fermionic states (trions, skyrmions,...) may be formed 
in a QHF under certain circumstances. 11-18 Magnetic impurities are characterized by their 
own spectrum of local spin excitations, and one can anticipate a strong interplay between 
local and collective excitations in a magnetically doped QHF. 

It is known that the influence of impurities on the discrete spectrum of quantized Landau 
electrons in a 2DEG has many specific features. Even such a basic property, as the interaction 
of a 2DEG with neutral short range impurities is far from being trivial. 19-21 Only the Landau 
states with a finite probability density on the scatterer locations interact with impurities. 
This means that the whole set of Landau states breaks down into two groups: the major part 
of the Landau levels (LLs) is not affected by the impurity scattering, and the states having 
a nonzero scattering amplitude on an impurity form a separate system of bound Landau 
states in the energy gaps between the free LLs. 

To be more specific, we consider a 2DEG formed in the n-type GaAs/GaAlAs heterostruc- 
tures doped with transition metal (TM) impurities. The reason for this choice is that the 
technology of (Ga,Mn)As epilayers is well developed, and the QHF regime is achieved ex- 
perimentally in GaAs based heterostrures. As a rule, transition metal ions substitute for 
the metallic component of the binary II- VI and III-V semiconductors. 22-24 The influence of 
isolated TM impurities on the spectrum of the Landau states was investigated in Ref. 25. It 
was shown that the resonance scattering in the ci-channel is in many respects similar to that 
of the short range impurity scattering in the s-channel. 19-21 The symmetry selection rules for 
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the resonance ci-waves in the cylindrical (symmetric) gauge pick up the Landau states with 
the orbital number m = (in the symmetric gauge). These states are the same states that 
are involved in the s-scattering by the impurities with a short range scattering potential. 21 
Besides, this scattering is spin selective in magnetically quantized 2DEG. 

It should be emphasized that in the problem under consideration the criterion of isolated 
impurities acquires a specific feature. In fact the Mn concentration range, where our theory is 
applicable, is limited from below by technological capabilities and from above by the obvious 
requirement that the impurity induced disorder does not destroy collective excitonic states. 
So, the relevant interval of bulk Mn concentrations is 10 13 cm -3 < n Mn < 10 15 cm -3 . Here the 
upper limit corresponds to the 2D concentration of 10 9 cm" 2 which in our case is actually well 
below the Landau band capacity at B ~ 10T that equals to the the 2D electron number 
on the upper (half-filled) LL. One may expect that the above mentioned classification of 
excitonic states is valid only at the bulk concentration n Mn < 10 15 cm~ 3 , which is much less 
than in the materials used for creation of dilute magnetic semiconductors. 26 

We calculate in this paper spectra of bound and continuous collective excitations related 
to magnetic impurities. When studying the influence of magnetic impurities on the excitonic 
spectrum of 2DEG, a distinction between the negative and positive signs of the gyromagnetic 
ratio of 2DEG electrons g 2DEG should be also mentioned. It will be shown that in the 
conventional situation of negative g 2 DEc the interaction with magnetic impurity lowers the 
ground state energy due to effectively antiferromagnetic character of the effective indirect 
exchange. This results in formation of a set of bound and delocalized collective excitations 
presenting combined modes classified by a change in the total spin number S z . When g 2 DEc > 
0, so that the g factors of both subsystems (2DEG electrons and impurities) have the same 
sign, magnetic impurities may form bound states in the gap below the spin exciton continuum 
and even initiate a global reconstruction of the QHF ground state. 

II. MODEL HAMILTONIAN 

Following Ref. 23, we describe the electron scattering on a TM impurity in semiconductor 
within the framework of the Anderson impurity model Hamiltonian 27 generalized for the case 
of multicharged impurity states in semiconductors. 28-30 According to this model, the principal 
source of magnetic interaction is the resonance scattering of conduction electrons on the d- 
electron levels of TM impurity in the presence of a strong on-site Coulomb interaction U. Due 
to this interaction, the local moment of TM impurity survives in the crystalline environment, 
and 'kinematic' indirect exchange interaction between the conduction and impurity electrons 
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arises in the second order in the s-<i-hybridization parameter, even in the absence of a direct 
exchange. 

The generic Hamiltonian describing the QHF regime in a magnetically doped semicon- 
ductor is 

H = H d + H s + H t + H sd . (2.1) 

Here H d = £V H di describes the TM impurities on the sites i, H s is related to the band elec- 
trons on the LLs, and H t is responsible for hybridization between the impurity ^-electrons 
and Landau electrons. Eventually, it is this hybridization that generates coupling between 
collective modes in 2DEG and localized spin excitations on the impurity sites. In our ex- 
tremely weak doping regime both the direct and indirect interactions between magnetic 
impurities are negligible. Each magnetic scatterer may be considered independently, and it 
is convenient to choose the symmetric gauge A = (—■§?/, \ x, 0) with the quantum numbers 
A = (n, m) for the Landau electrons hybridized with the atomic ci-electrons centered around 
the site i positioned in the center of coordinates. The Coulomb interaction is taken into 
account in the impurity and in the band electron subsystems. Besides the direct Coulomb 
interaction between the d and s electrons described by the last term in the Hamiltonian (2.1) 
is added to the conventional impurity Hamiltonian (cf. Ref. 27) described by the first and 
third terms. All additional interactions will be discussed below in detail. 

Substitutional Mn impurity in GaAs retains all its five d electrons due to a special sta- 
bility of the half filled 3d shell. In the p-type GaAs the electrically neutral state of Mn 
in Ga position is Mn( 3+ )(3<i 5 -|-hole), where the hole is bound on the relatively shallow ac- 
ceptor level near the top of the valence band, whereas the occupied ci-electron levels are 
deep in the valence band. 23 ' 24 ' 31 ' 32 In the n-type heterostructures these acceptor states are 
overcompensated, and the chemical potential is pinned on one of the lowest Landau levels in 
the conduction band. Since we are interested only in the low-energy excitations above the 
ground state of n-type system, Mn impurities will be considered as the Mn( 3+ )(3<i 5 ) ions in 
the subsequent calculations. 

A. Single-orbital model. Spin-selective hybridization. 

One may significantly simplify the calculation of the spectrum of excitations by reducing 
the general Hamiltonian (2.1) to the form, in which only the terms relevant to the calculation 
of desired collective states are present. As a result of this simplification outlined below in 
Subsection II B one arrives at the single-orbital, single Landau band hybridization Hamilto- 
nian, which explicitly takes into account the Hund rule governing the high-spin states 3d 5 
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and 3c? 6 of the Mn 3d shell (the state 3c? 4 is proved to be irrelevant in our specific case of 
Mn in GaAs lattice, see below Fig. 1). These impurity states are characterized by the max- 
imum total spin quantum numbers = 5/2 at 3c? 5 , and = 2 at 3c? 6 , and the effective 
Hamiltonian H is defined in the charge sector {|c? 5 ,iV), \d 6 ,N — 1)} of states with variable 
number N or N — 1 of the electrons on the highest n-th LL (of course in our case N ~ N^). 
The Hamiltonian now reads 

H = e dCT n 7()CT + U n 7oT n 7oi + ^ e ntJ at nmc a nrna + H' Coul + H t . (2.2) 

a ma 

where the impurity Hamiltonian of Eq. (2.1) is represented in by the two first terms, 
in which n 7()Cr = c 70O .c 70(T , and c 7()Cr is the creation operator for the o?-electron at the orbital 
7o with the spin z-component a. The notation 7 designates the only 3c?-orbital with the 
Y 02 ~ 3z 2 — r 2 symmetry, which effectively couples with the m = state of the LL. 25 The 
parameter U characterizes Coulomb and exchange interactions determining the addition 
energy for the transition 3c? 5 — > 3c? 6 . The third term in Eq. (2.2) is the Hamiltonian 
of noninteracting Landau electrons where (v nma is the creation operator for the (n, m, a) 
Landau state. The most of interaction components are included in H' Cou i- This term does 
not include only the d-d interaction parametrized by U and the last term H t . The latter 
generically is also the part of Coulomb interaction between impurity and Landau electrons 
which intermixes impurity and Landau orbitals . However, in our case H t acquires the form 
of single-electron hybridization operator [see discussion after Eq. (2.13)], 

H t = J2 W noal 0a c Joa + H.c. (2.3) 

a 

As was mentioned above, this operator is responsible for the resonance orbital-selective scat- 
terings in QHF. It includes hybridization of the impurity electron with the 2DEG electrons 
within the n-th LL. This means that only the influence of impurity on the intra-LL excita- 
tions (of the spin- wave type) is taken into account. The hybridization with the states with 
n' 7^ n describing the processes with energy change hw c or higher is omitted. 

In the absence of interaction term H' Coul , the Hamiltonian (2.2) acts in the subspace 

|c? 5 , s; vac), |c? 6 , s + ^; a n0T |vac) and |c? 5 , s + 1; 0^0^! vac), (2.4) 

where the fully polarized 2DEG without impurity is chosen to be the 'vacuum' state |vac) =| 
I). Therefore aj^lvac) = a ram j|vac) = 0. We represent the total spin component as 
S z = ~y + s. Thus we characterize the states in the set (2.4) by the quantum number S z . 
It is important that only S z = S z ^ + is an exact spin quantum number in our system. 
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Separately, the Hamiltonian (2.2) commutes neither with the spin component of the LL 
electrons nor with the impurity spin component Si d \ Equally, it does not commute with 
the total spin S 2 and with the spins (S 1 ^) 2 and (S^) 2 (see Appendix B). The number s 
in the set (2.4) changes within the interval — | < s < |. It is convenient to choose the 
state |c/ 5 ,|;vac) as a reference point ('global vacuum'). This state is not mixed with any 
other state of the system by the operator (2.3) so that it enters the set of eigenstates of the 
Hamiltonian (2.2), although at #2Deg < it is one of the excited states of the system. 

Within a given 'triad' (2.4), i.e. at a given s, the operator (2.3) intermixes these basis 
states. The corresponding non-diagonal matrix elements are (vac; s, d 5 \H t \d 6 , s + |; a n0 tl vac ) 
and (vacla^; s + \, d 6 \H t \d 5 , s + 1; a* , a n of|vac), where the bra- and ket-vectors are appro- 
priately normalized. Therefore, for any given quantum number S z = ^ + s the magnetic 
impurity scattering problem can be effectively described in terms of a single- orbital impurity 
model that involves only one or two G?7 -electrons. The single-orbital basis 

|s_;vac), \s ; a n0T |vac) and \s + ; al ol a n0] \vac) (2.5) 

arises instead of the original multi-electron basis (2.4) where the indices (— , 0, +) label the 
bare energies E s _, E So = E s _ + U + e dT - £„ T and E s+ = E s _ + (^ - g 2 nv G )HBB. Here 
giUsB = e^t — e di, and ^deg^b-B = £ n | — £ ni are Zeeman energies for impurity and 2DEG, 
respectively. The two states E s± form a Zeeman doublet for a given s, while the state E s0 
becomes resonant with the LL continuum. 

As a result of this mapping, where only one component 7 of the normalized multi-electron 
states is responsible for the hybridization, the H t operator can be redefined for each triad 
(2.5) as 

H t (S z ) = V [/3 T ( S )c{a + /3i(s)c{&o] + H.c. (2.6) 

It becomes thereby spin-selective. The shorthand notation ao = a n0 i, bo = a„ | is used here 
and below; cj^ is the creation Fermi operator for the s± impurity states \s±] vac) = cjyjvac) 
and |s ;vac) = c|c||vac). The Clebsch-Gordan coefficients f3 a {s) reflect normalization of 
eigenvectors (2.4) by replacing them with normalized single-orbital basis (2.5). For s = 
-5/2, -3/2, -1/2, 1/2, 3/2 we have 



The highest state in the bare Zeeman ladder | (|)_;vac) = \d 5 , |;vac) remains nonhy- 
bridized. 

Unlike the original Anderson model 27 the mixing coefficient, V = W n o in our particular 
case, arises as a non-diagonal component of the s-d Coulomb interaction (see Subsections 
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IIC and II D for further discussion). 

B. Description of the employed simplifications 

Here we list the simplifications which have allowed us to reduce Eq. (2.1) to the Hamil- 
tonian (2.2) with following change H t — > Tt t , and to apply to our system. 

The first simplification exploits the fact that the characteristic Coulomb energy of Landau 
electrons Eq = ae 2 / kIb is considered to be small in the QHF regime as compared to the 
cyclotron energy Jtw c . Here a is the average form-factor related to the finite thickness of 
the 2DEG (0.3<a<l). In the E c <ti fku c limit one may neglect the LL mixing. Besides, it 
implies that in our case the energies of collective excitations are smaller than fuv c . 




(a) (b) 

FIG. 1: Allen reactions that involve an additional electron (a) or hole (b) in the impurity site. 
edi and + U are addition energies for the 5-th and 6-th electrons in the 3d-shell of Mn ion in 
accordance with Eqs. (2.8). The ground state with the energy Eq corresponds to the completely oc- 
cupied lowest Landau subband. Spins of Mn 3d shell and occupied Landau subband are antiparallel 
because of the different signs of g factors for Mn and 2DEG in GaAs. 

The second simplification is related to the 'deepness' of the 3d- levels of a neutral Mn 
impurity relatively to the bottom of conduction band in GaAs. We know from the previous 
studies 25 that the scattering potential created by a TM impurity for the Landau electrons 
is generated by the s-d hybridization. It has the resonance character, and the spin selective 
scattering becomes strong when one of the impurity 3d levels is close to the LLs of conduction 
electrons. The process of s-d hybridization may be represented by the so called 'Allen 
reactions' 23 ' 24,33 (see Fig. 1) 

3d 5 ^3d 6 + h, (2.8a) 
3d 5 ^3d 4 + e. (2.8b) 



8 



The first of these reactions describes hopping of an electron from the filled Landau subband 
to the impurity rf-shell, whereas the second one means hopping of an electron from the d 
shell to a state in the empty Landau subband. It is known from the numerical calculations 32 
that the addition energy for the 6-th electron in the Mn 3d shell (e^ FR state in terms of 
Ref. 24) is in resonance with the states near the bottom of GaAs conduction band. It really 
means that the values U and e d ^ well compensate each other in the sum + U. So, one may 
retain only the processes (2.8a) in H t and neglect contributions of the 3d 5 — >■ 3d 4 ionization. 

The third major reduction of the Hamiltonian is the elimination of the impurity orbital 
degrees of freedom due to the selection rules for the s-d hybridization matrix elements. 25 This 
orbital selectivity arises, first, because of symmetry reasons since only electrons with equal 
axial m-numbers in d and LL states are hybridized. Second, a precondition of the selectivity 
is related to the fact that the magnetic length l B is much larger than the radius r d of 3d- 
electron state (in the energy scale this condition takes the form of inequality U, e da fajJ c )- 
The hybridization integral determined by the overlap of the d- and Landau wave functions 
behaves as ~ (r d /l B ) m C 1 for m ^ 0. All resonance scattering (hybridization) amplitudes 
with m 7^ are thus negligibly small, and only the s-scattering term (m = 0) can be 
retained in H t . This explains why only one of the five 3d-orbitals, namely 70, is involved in 
the resonance interaction with the 2D Landau electrons. 



C. Interaction Hamiltonian in excitonic representation 

As it was mentioned above, the states of the system are characterized by the total spin 
component S z . For a given S z we may deduce the effective Hamiltonian 

H{S Z ) =H d + H[ s) + H t {S z ) + H„ + it-d (2.9) 

with the single-orbital impurity term H d = e^t^t + e d|^-| + Unin± (n a = cJ.c CT ), and with the 
hybridization term determined by Eq. (2.6). 

The remaining terms in the Hamiltonian (2.9) are defined within the framework of the 
single-LL approximation for the Landau electrons. 3_6 ' 9_12 ' 16 ' 18 ' 34-36 Although only the states 
with m = in the LL are involved in the resonance scattering the complete basis for the 
description of collective excitations includes all m-orbitals of the LL, and the corresponding 
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Schrodinger field operators should be taken in the form 

m 

^(R) = c^ d (R) + C.(2)X]^m(r). (2.10) 

m 

Here the shorthand notation a m = a nm ^ b m = a nm i is used. R = (r, z) is the 3D coordinate 
with the reference point at the impurity site, ( s {z) is the size-quantized functions of s- 
electrons in the layer, ip m is the wave function of the n-th LL, where index m in the symmetric 
gauge changes within the interval (—n, —n + 1, ...N^ — n — 1). 

Using the above definitions and Eqs. (2.10) in the generic interaction operator 

^coui = \ J d 3 R 1 d 3 R 2 ^l 2 (R 2 )^ ai (R 1 )W(R 1 - R 2 ) * (Tl (R 1 )i CT2 (R 2 ) (2.11) 

[where W(R) ~ e 2 / kR at R ^> r d ], one may rewrite the s-s and s-d Coulomb interactions in 
the excitonic representation (ER). 6 ' 9 ' 18 This actually means that after substitution of Eqs. 
(2.10) into formula (2.11) the latter can be expressed in terms of combinations of various 
components of the density-matrix operators. These are so-called ER operators presented in 
our case only by the intra-LL set, i.e. by the spin-exciton operators <2 q where an electron 
is promoted from one spin-sublevel to another [see Refs. 6,9,18 and Appendix A where the 
necessary ER equations are given with the reference to our case] and by operators A^ and 
acting within the sublevels a or b (see ibidem). As a result the Coulomb terms of Eq. (2.9) 
can be written only by means of the intra-sublevel operators A^ and £> q [their definitions 
are given by Eq. (A. 4) in Appendix A], 

= ^ E w »b) ( A i A « + 2A i B « + B i B *) - \ + Bo) J2 w s-M) , (2.i2) 

q q 

= (n T + fii) J2 Wsdiq)^ + B q ) . (2.13) 
q 

The Coulomb vertices are presented also in Appendix A [Eqs. (A12) and (A13)]. 

We neglect in Eq. (2.11) the direct exchange s-d interaction terms (see the next sub- 
section). The mixing operator H t (S z ) in our model Hamiltonian (2.6) includes in fact off- 
diagonal interaction terms from Hcoui- Indeed, Coulomb interaction described by the terms 
~ cjn-f& m + H.c. and ~ cjnj<2 m + H.c. induces transitions adding or removing one electron 
to the d-centei in accordance with the Allen reaction diagrams (2.8). These terms represent 
the s-d hybridization formally conditioned by the ci-center occupation; however, since in our 
case the reaction (2.8b) is forbidden, they actually operate as ~ c\b m +H.c. and ~ cla m 
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+ H.c. in Eq. (2.6), respectively. [In terms of the d 5 <-> d 6 transitions the hybridization is 
taken just in the form of Eq. (2.3).] 

The single particle Hamiltonian for LL electrons may be also written in the ER represen- 
tation, 

H{ s) = {(e n - e z /2)Ao + (e„ + s z /2)B ] , (2.14) 
where e z = I^degI/Ab-B and e n = (n + l/2)hu c . 

D. Numerical estimates of the energy parameters 

Before turning to our main task, i.e. to the calculation of excitation spectra, it is worth- 
while to evaluate the characteristic energy parameters related to this problem. We estimate 
the parameters of 2DEG in GaAs for the typical value B = 10 T of magnetic field. In 
this field Eq ~ 5meV characterizes the Coulomb interaction (A13). Below in our calcu- 
lation this value is mostly presented by the spin-exciton mass, which can be estimated 
empirically, i.e. the inverse mass is 1/M X ~ 2meV in energy units. The LLs' spacing is 
fkv c ~ 16meV, and the Zeeman splitting between two Landau subbands is e z ~ 0.25 meV 
(because g 2DEG = g GaAs ~ —0.44). The Zeeman splitting for Mn ion is Qi^bB ss 1.1 meV 
(because gi = g Mn ~ 2.0). The hybridization constant V and the repulsion U are the other 
important parameters characterizing the magnetic impurity. It is rather difficult to extract 
them from the available experimental data. We can only roughly estimate the energy U 
as a distance between the Mn-related peaks in the density of states of occupied and empty 
states in the spectrum of bulk (Ga,Mn)As calculated with an account of the electron-electron 
interaction. 32 Such an estimate gives U ~ 4 — 4.5 eV. From the same calculations we estimate 
the energy difference 

A = e d ^ + U-e n + e z /2, (2.15) 

which determines the position of the Mn(<i 6 ) electron level above the bottom of Landau 
band (see Fig. la) as A < 0.1 eV. In order to estimate the parameter V, one should recollect 
that the dominating contribution to hybridization integral is given by the matrix elements 
of Coulomb interaction, having the form Vcjaocjcj (see above). This means that V ~ 
Ur^ 2 ((z d )/l B . Estimating the radius of the ipj function as r d ~ 2 A, and ((z d ) ~ 0.15 A -1 / 2 
(for the impurity located in the vicinity of the quantum well bottom), one gets V ~ 20 meV. 
This gives | V| 2 /A ~ 4 — 8 meV for the relevant kinematic exchange parameter. At the same 
time the direct exchange turns out to be insignificant. Indeed, one can estimate from Eq. 
(2.11) that the characteristic coupling constants for the terms ~ ctcjfe^a^ and ctc\b mi a m2 
are of the order of Ur d \((z d )\ 2 /l^ being therefore by the factor ~ A/U smaller than |V| 2 /A. 
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III. COLLECTIVE SPIN-FLIP STATES. NEGATIVE gf 2 DEG-FACTOR 

The Coulomb interactions H s - S and H s ^d admix the LL states with m ^ to the three- 
state basis (2.5). Instead of triads (2.5), the basis 

|s_;vac), |so;a-m|vac) and \s + ; <2 q |vac) (3.1) 

contains spin-exciton continua <2 q |vac) attached to the spin-flipped impurity state s + . [The 
definition of the spin-exciton creation operator is given by Eq. (Al).] 

This set is complete only within the single-orbital approximation. 37 At a given s it is 
convenient to take the energy E 0+ (s) of the state |0) = |s + ;vac) = c||vac) as the reference 
point because this state is not affected by the hybridization within the framework of the 
single-orbital model. This energy is defined as E 0+ (s) = (vac|;s + 1, d 5 \H\d 5 , s + 1; |vac) 
where the Hamiltonian H is given by Eqs. (2.2) and (2.3). For a given S z — ^ + s we have 
E 0+ (s) = E vac — (| — sjgiHsB with E vac defined as the energy of the global vacuum state 
\d 5 ,5/2; |vac). One can check with the help of expressions (Bl) and (B4) in Appendix B 
that the vectors \sq; a m |vac) and Q q |0) at q ^ correspond to the definite total spin state 
with S = S z , whereas |s_; |vac) and Q o \0) are not characterized by any definite number S*. 38 



A. Secular equation 

Following the above discussion the spin-flip operator may be represented in the form 

m q 

The normalizability condition (X\X) < oo for the bound spin- exciton state \X) = X*\Q) 
then reads J2 m \D m \ 2 + X] q Ifi^l 2 < 00 an d the sum 

A fe = El/(q)l 2 = §/^l/(q)l 2 

q 

presenting the contribution of continuous spin excitons into the norm (X\X) becomes thereby 
an essential characteristic of the spin-flip excitation. For the regular (normalizable) solutions 
we expect /(q) ~ A^ 1 ^ 2 . Besides, singular states, for which the sum J^cilfi^l 2 diverges 
also exist. These states form continuous spectrum of impurity- related spin-excitons. 
The coefficients D m and /(q) are determined from the equation 

[H,Xi]\0)=E\X). (3.3) 

where the energy E is counted from E 0+ (s). Before turning to the computation we specify 
the energy levels of the basis states (3.1) at V — 0. The state |s_;vac) has the energy 
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E _(s) = E 0+ (s) — gifi B B. The doubly occupied impurity state \d ; a m |vac) appears due to 
the charge transfer with creation of a conventional 'hole' in the LL. Its energy is E d:tn (s) = 
Eo+ + £d, m where 

£d,m = e dl + U + e z /2 - e n + e m + (3.4) 

[cf. Eq. (2.15)]. Here = (l/JV^qW,,^) [see Eq. (All) for definition of £«,]. This 
term appears due to the global electroneutrality requirement when calculating the energy 
of the hole a m |vac). 4 ' 5 ' 9,11 The term e m = -(2/N^,) £ q h m+rit7n+n (q)W sd (q) is the Coulomb 
interaction energy of the hole a m |vac) with the doubly occupied <i-center [see Eq. (A2) for 
functions hit}- 

Substituting operators (2.9) and (3.2) into Eq. (3.3), projecting the result onto the basis 
vectors (3.1) and using Eqs. (2.6), (A3), (2.12), (2.13), (2.14) and (A5)-(A10) we obtain a 
closed system of equations for the coefficients D m , /(q). This system defines the eigenvalues 
of Eq. (3.3) for a given s. The symmetry of the problem allows us to look for the solutions 
in the form /(q) = f m (q)e im<1 '. Below we limit ourselves to a study of the isotropic case 
of m — 0. (Discussion of the case m ^ may be found in Ref. 35.) As a result we get 
D m = D S m fi and our system for given S z = ^ + s acquires the simple form 

E + 9il i B B = Pi(s)V*D , 

(E-£ dfi )D = fy(s)V - /^(s) VJV~ 1/2 ^2 KniQ) f(v) . 

q 

(E-e z - S q ) f(q) = -N^ 1/2 h nn (q)^(s)V*D . (3.5) 

The energy of the free exciton state Q q |0) is e 2 + £ q [see Eq. (All)]. 

The collective states localized around a magnetic impurity are described by solutions of 
Eq. (3.5) outside the free spin-wave band (i.e. in the energy interval E < e z or E > £z+£oo)- 
The corresponding eigenfunctions are characterized by the regular envelope function fo(q). 
We arrive then at the secular equation 

g(f) v Kn(q)\ 2 W _ E-S dtQ 

Nej, e — sz — Sq E + g lf i B B~ \V\ 2 l ' J 

for the energy E. The first term in the l.h.s. of Eq. (3.6), including the sum of spin-exciton 
propagators, presents the self energy, which usually arises in the Schrodinger or Lippmann- 
Schwinger equation describing the perturbation introduced by a short-range potential into 
the continuous spectrum. The prototype of this term in the theory of magnetic defects is 
the self energy for localized spin waves in the Heisenberg ferromagnet with a single substi- 
tution impurity. 39 Specific features of our model are manifested by the energy dependence 
of impurity related processes. First, instead of a constant term (inverse impurity potential) 
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in the r.h.s. of Eq. (3.6) we have the inverse resonance potential \V\ 2 /{E — £4,0), which 
stems from the hybridization between LLs and the 3d-level of impurity electron. 28,29 Second, 
an additional term describing impurity spin-flip process in terms of the single-orbital model 
arises in the l.h.s. of Eq. (3.6). 

B. Spectrum of the localized states 

First we carry out a simple study considering solutions of Eq. (3.6) in the absence of an 
exciton band, i.e. by formally substituting £ q = into Eq. (3.6). (This is instructive in 
order to classify the bound collective states.) We obtain then a simple algebraic equation 
with two roots. When solving this equation we use the sum rule \h nn (q)\ 2 = and 
neglect the energy dependence in the r.h.s. due to the condition £ dj0 ?a A 3> E. Each 
doublet is bound to its own reference energy Eq + (s) in accordance with the corresponding 
spin component = s + 1 of the Mn' +2 ' ion. Due to the kinematic exchange (2-nd order 
spin-flip processes) each state in the Zeeman grid (lower root of Eq. 3.6) acquires a partner 
state (upper root of Eq. 3.6), except for the highest level with s = 5/2, which remains 
intact, because the spin flip processes are kinematically forbidden for this state. The level 
splitting is illustrated by the scheme in Fig. 2. We see that the kinematic exchange makes 
the Zeeman states of impurity ion non-equidistant and an additional multiplet of excited 
states arises as a prototype of the bound spin-excitons. 

Having this classification in mind we turn to calculating the bound exciton states for a 
finite dispersion of the free spin waves. According to the estimates of the model parameters 
presented in Subsection II C we solve Eq. (3.6) for the realistic conditions E c > 9%^bE 3> £z 
whereas the ratio between the energies E c and |V| 2 /A may be arbitrary. 

All the generic features of impurity-related states may be seen in the case of a unit filling 
where n = (is — 1) and we study this situation in detail. The solutions we are looking for 
are localized in the energy interval \E — £2,0 1 ~ A so one can neglect the energy dependence 
in the r.h.s. of Eq. (3.6). A graphical solution of Eq. (3.6) is schematically shown in Fig 3. 

Two intersection points labeled as E^ and E^ correspond to two discrete solutions. Just 
as in Fig. 2 this pair of solutions arises at any s except for s = 5/2. The lower solution with 
the energy E^ is the state of the Mn' +2 ^ ion with the spin component (S^) ~ s shifted 
downwards from the value E vac + (s — ^)gifj,sB by the effective exchange interaction with 
the spin- wave continuum (in this case D > iV b ). The upper solution corresponds to the 
spin-flipped state of the Mn' +2 ' ion with (S^) « s + 1 dressed with the spin- wave localized 
on the impurity. In this case (S^) ^ — 1 and N b > D . This bound exciton state, 
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FIG. 2: A scheme of the Mn' +2 ' Zeeman level splitting due to the kinematic exchange in the absence 
of exciton dispersion. The bare Zeeman ladder is shown on the left. Five of six levels in this grid are 
shifted down (extreme right column), whereas the s = 5/2 level remains not renormalized. Each of 
these five levels has its high-energy counterpart. The energy is measured in the gi^sB units. The 
following values of the input parameters are chosen: ez = 0.2, |F| 2 /A = 2 and £d,o = A>1. The 
factors Pf(s) and /3 2 (s) are presented by Eq. (2.7). 

described semi-phenomenologically in Ref. 35, is shallow compared with the main charac- 
teristic energy parameter E c . Like in many other impurity- related states in 2DEG 21 ' 25 ' 40 its 

(s) 

energy is confined within the interval —gifisB < < £z in the logarithmic vicinity of the 
bottom of the delocalized spin-exciton band. Due to this fact one may find the level position 
analytically. Using the quadratic approximation for the exciton dispersion law £ q = q 2 /2M x 
and turning from summation to integration in the l.h.s. of Eq. (3.6) one has 

^ 2 v- IM<?)I 2 



/3fM x ln[ 7 M x (|ez-£|)] (3.7) 



Here M x is the spin-exciton mass denned as 1/M X = J °° dpp 3 v ss (p)e p2 / 2 /2 ~ E c [see Eqs. 
(All) and (A13), the = ! units are used.] and 7 = 1.781.... Then the binding energy 



2 



£ *" = £z " w, exp (-#i,B " W^n< J ' (3 ' 8) 

is found from Eq. (3.7). This result is valid provided at least one of the two inequalities, 
(3^M x gifi B B <C f3 2 or (3 2 M X \V\ 2 / A <C 1, holds, which is not too strict requirement due to 
the exponential smallness of the second term in r.h.s. of Eq. (3.8). 

The asymptotic value of the lower state is also easily found. In the case of strong 
hybridization |V| 2 /A ^> Ec one gets E^ —QifisB — 5|y| 2 /6A. In this asymptotic limit 
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FIG. 3: Graphical solution of the secular equation. The l.h.s. and r.h.s. of Eq. (3.6) are shown as 
functions of the argument E by solid and dashed lines, respectively (-E-dependence in the r.h.s. is 
neglected). Filled area indicates possible values of the l.h.s. because it belongs to the interval of E 
where the sum in Eq. (3.6) becomes indefinite. 

the excitation energy does not depend on s. In the opposite limit | V"| 2 /A <C E c we have 
E ( S) ^-g^ B B-(3 2 {s)\V\ 2 /A. 

In the intermediate region \V\ 2 /A ~ E c Eq. (3.6) for E { s) can be solved numerically. It 
is convenient to rewrite this equation in the dimensionless form 

fl8) fi^. + |!) +1 = o (3.9) 

where £ = A/M X |T/| 2 is the ratio of the characteristic Coulomb energy in the Landau band 
and the characteristic kinematic exchange energy. The relevant energy parameters in (3.6) 
are redefined as E = (\V\ 2 /A) F^(^), e z + £ q = M~ x e{q) and g,^ B B = (\V\ 2 /A) g. Then 

~ (s) 

the system of localized levels Eq£ counted from the global vacuum energy is described by 
the set of equations 

E ( o,l = -giUBBW -s) + F$(0\V\ 2 /A (3.10) 

with s = -5/2, -3/2, -1/2, 1/2,3/2. The family of lower roots F (s) (£) of Eq. (3.9) changing 
smoothly from —j3 2 — g at £ = oo to approximately —6/5 — 5g(3 2 /6 at £ = describe the 
renormalization of the Zeeman grid of impurity spin-flipped states due to the kinematic 
exchange with LL continuum. To illustrate this dependence we have found the solution 
of Eq. (3.10) for g = 0.25 neglecting e z and modeling the spin-exciton dispersion by the 
function e(q) = 2 — 2e~ q2 ^I (q 2 /4), which corresponds to the ideal 2D case. 3-5 (At the 
same time the parameter M x may be considered as an empirical value.) The results of this 
calculation are presented in Fig. 4. 
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FIG. 4: The lower root of Eq. (3.9) with g = 0.25. The numbers s are indicated near the curves. 
See text for further details. 

C. Delocalized impurity-related excitations. 



We conclude this section by a brief discussion of the delocalized states (free spin waves 
distorted by the resonance magnetic impurity scattering). These states are described by the 
functions f(q) with a divergent norm in the expansion (3.2). Secular equation for these states 
cannot be presented in the form (3.6) but there are solutions satisfying Eqs. (3.5) at any 
energy within the spin-exciton band, Ez < E < ez + £oo- These states are the 'counterparts' 
of the levels E^ in the spin wave continuum. Let qo(E) be a root of equation ez + S qo = E. 



Substituting 



(3.11) 



into Eqs. (3.5) one gets three equations for the coefficients Do, C and u(q). Turning 
from summation to integration and using the rule X) q ^|q|,go = ^qoL/n (L 2 = 2izN ( f > , being 
the 2DEG area) one finds the coefficient u(q) from the equation (3 2 (E — ez — £ q )u(q) = 
-(3lN~ 1/2 h nn {q){gifi B B + E). Then equation 



C e -?o 2 /4 = 1 + 



(E + gifi B B) 







f qdq\h nn \ 
J E-e z - 



\ 



+ 



S. 



d,0 



E 



V\ 



\ 



(3.12) 



/ 



for the spectrum is derived from Eq. (3.5) in the thermodynamic limit (L, — > oo). It can 
be readily seen that the norm of the function (3.11) diverges as ^2 \u(q)\ 2 ~ N^. 
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IV. POSITIVE gf 2DEG -FACTOR. PINNING OF THE QHF SPIN 

Experimentally the magnitude of the g 2DEG factor in GaAs/ALGai-^As structures can be 
altered gradually by changing pressure or by varying Al content (x). It can be made very 
small and even change its sign. 14 The value of # 2 deg m ay be effectively reduced also due 
to optical orientation of nuclear spins changing the electron Zeeman splitting (Overhauser 
shift). 13,17 In this section we discuss the impurity-related reconstruction of the ground state 
and the spectrum of spin-flip excitations at small but positive values of g 2 DEo- It will be 
shown below that even a minute amount of magnetic impurities can drastically influence the 
QHF state. 

Keeping the previous notations, it is now convenient to redirect the magnetization axis 
(z — > — z), i.e. to make formal transformation ^ — > — gi instead of changing the sign of 
(?2deg- It is clear that at least in the absence of the s-d hybridization the global vacuum 
state \d 5 , 5/2; |vac) serves as the ground state, and all the spin-flips cost positive energy. 
The localized states can still be found from Eq. (3.6) with redefined Zeeman energies, 
giUsB — > —gifi B B and e z — > e^. The latter parameter actually takes the values e*^ = 
92DEgI jL bB>0.1K. Making change g — > — g in Eq. (3.9), we denote the lower root of this 
new equation as F^- This root corresponds to the energy of the localized spin exciton with 
changed impurity spin projection, 8S^ ~3/2 — s, where s= |, |, — |, — |, — |. The total spin 
component is S z — — — s (when presenting results, we return to the 'normal' coordinate 
system where z is directed along B), and the energy counted off the global vacuum level is 
given by formula 



(It should be noted that now the new global vacuum is really below the old one by the energy 
bgiHsB.) Functions F^\C) are presented in Fig. 5. 

Other roots of the secular equation belong to the continuous spectrum. These states may 
be analyzed following the approach described in Subsection III C. The special 'resonance' 
solution of Eq. (3.12) with substituted for — gi arises in this case at E = gifi B B > e* z . Then 
u(q) = and the norm (X\X) diverges not as ~ but as L ~ N^ 2 (see discussion in the 
next section). As a function of s, the delocalized 'resonance' states form a set of equidistant 
levels 



(again the energy of global vacuum is taken as the reference level). 

When looking for the x-type solutions at E <e1 L but \E\ 1/M X , one may use Eq. 
(3.7). Then one obtains for the localized spin- exciton energy 41 E^ = |\/| 2 Fi s) (£)/A the 



= g^ B B(3/2 - s) + i^(0|Vf/A . 



(4.1) 




(4.2) 
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FIG. 5: The lower root of Eq. (3.9) with negative parameter g. Calculation is performed for 
g = — 0.25 and e(q) = 2 — 2e~ q2 ^I^q 2 / A). The values of s are indicated near the curves. 



following equation 



7 M X 6XP \pfM x g tf i B B PfM x \V\< 



ft 



(4.3) 



instead of (3.8). Here s — § has to be taken for the excitation from the ground state, then 
fll — | and (3 i = l. The exponentially small energy E^^ thus corresponds to formation of 
bound spin exciton of large radius. However, for sufficiently small e\ (or for a strong enough 
kinematic exchange) the energy E^ 2 ^ becomes negative which means an instability of the 
global vacuum \d 5 ,5/2; |vac) considered as the QHF ground state. This instability appears 
provided 

£ < U , (4-4) 

where £ c i is determined by the equation 



\v\ 2 eW 2 \u)/a + s* z = o 



(4.5) 



The question arises, whether the condition (4.4) mean the global reconstruction of the 
ground state and appearance of a new state with many spin excitons bound to the magnetic 
impurity? To clarify this point, we discuss the limiting situation where e\ — > but still 
N^e^^oo. Then the ground state at any £ is no longer the global vacuum because creation 
of one spin exciton bound to the impurity lowers the energy of the system. The corresponding 
energy gain compared to the global vacuum is presented as \V\ 2 G 1 (^)/ A. [The subscript 
'1' corresponds to one bound exciton; specifically, we have Gi(£) = F^ 2 \] To answer 
the question, one should consider the situation with K captured spin excitons (then S z = 
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K ^— ) and calculate the proper value Gk{£) at arbitrary K. The latter is determined 

by the competition between the antiferromagnetic kinematic exchange, which forces 2DEG 
spins to reorient in the direction opposite to the impurity spin, and the Coulomb-exchange 
energy appearing due to 2DEG inhomogeneity in a cluster of K spin excitons bound to 
the impurity. This inhomogeneity energy is measured in 1/M X units. Calculation of Gk at 
K > 1 (but 1) is beyond the abilities of our present approach but we can consider the 
case of K ^> 1 and find the conditions under which such a massive pinning of 2DEG spins in 
the vicinity of the impurity turns out more advantageous than binding of single spin exciton 
(i.e. Goo>Gi). 



A. Skyrmionic states created by magnetic impurities 

The state with K >> 1 can be described as a collective topological defect (skyrmion) 
pinned to a magnetic impurity. 36 A smooth inhomogeneity in the system of spins may be 
presented as a continuous rotation in the 3D space. If one characterizes the local direction 
of the spins by a unit vector n(r) with components n x = sin 9 cos ip, n y = sin 9 sin <p, and 
n z = cos 9 ((p and 9 are the two first Eulerian angles) then the conditions 0| r=o = and 
0| r=oo = 7r inevitably result in the appearance of the topological invariant q T = f drp(r) 
where the density 

p(r) = — n ■ (d x n) x (d y n) (4.6) 

47T 

is a vortex characteristics of the spatial twist. The value q T has to be an integer nonzero 
number. 42 Its physical meaning is the number of excessive (q T < 0) or deficient (q T > 0) 
electrons in the system, 11 ' 12 ' 16 ' 18 i.e. q T = — N. In a perfect 2DEG and at nearly zero 
Zeeman gap (e% — > 0) such a weakly inhomogeneous skyrmion state has the energy 

3 1 

£sk = -£oo<?t + 2fif (I9t| - 9t) • (4.7) 

This result is valid within the single Landau level approximation (see, e.g., Ref. 18). It is 
enough to consider the case g T = ±l, because any state with |g T | > 1 is merely a combination 
of 'singly-charged' skyrmions. Due to the hybridization with the impurity the skyrmionic 
state gains a negative kinematic exchange energy. The latter has to be taken into account in 
combination with the Coulomb-exchange energy (4.7) and with the finite positive Zeeman 
energy at g* DEG > 

E z = e* z K, where K = J (1 + cos9)dr (4.8) 
(in the clean 2DEG the skyrmion energy is given by E^+Ez). 
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One can conclude from symmetry considerations that the impurity is located at the cen- 
ter of the topological defect. Then additional pinning energy may be found by means of 
the conventional energy minimization procedure where the Euler angles are used as varia- 
tional parameters. This energy is a difference between the energy of the global vacuum state 
with a distant skyrmion and the ground state energy calculated in the presence of magnetic 
impurity at the center of the topological defect (cf. Ref. 36 where similar procedure was 
elaborated in the limit of potential impurity scattering). Namely, to calculate the contribu- 
tion of magnetic impurity at K^$> 1, one should consider a domain around impurity which is 
small in comparison with a characteristic area of the skyrmion, but contains a large enough 
number of spin-flipped LL electrons involved in the formation of pinned topological defect. 
Then the situation becomes similar to that considered in Sec. Ill: s-d hybridization of the 
impurity electron with the m = electron in this domain generates the kinematic exchange 
in accordance with Fig. 1, and leads to reconstruction of the spectrum in accordance with 
Eq. (3.6). The shift of the energy with respect to the global vacuum is given by the value 
5gi/j,BB + EQ~ 5 ^ 2 \ where Eq~ 5 ^ 2 ^ is determined by Eq. (3. 10). 43 Hence we obtain that the 
pining energy is E sKpin = -gi/ji B B - F { ~ 5/2) (£) | V\ 2 /A, where F ( ~ 5/2) is shown in Fig. 4. In 
the limit of strong pining (£ , S k, p in £sk) and 'frozen' impurity spin (g ^> 1) this result agrees 
with the pinning energy found earlier. 36 

The energy -E s k, P in is calculated in the leading approximation, which does not depend on 
the charge q T . However, it is instructive to obtain the correction related to the inhomogeneity 
of the texture. It is known 16 ' 18 that the density (4.6) may be interpreted in terms of effective 
magnetic field appearing in the Schrodinger equation due to this inhomogeneity. One may 
introduce the renormalized magnetic length Ib ^ Ib as 

l = l-27rp(r). (4.9) 

l B L B 

Taking into account that |V| 2 /A ~ l/l B and £ ~ l B and rewriting Eqs. (4.6) and (4.9) 
in terms of the Euler angles, 42 one finds the correction to pinning energy due to the finite 
radius R* of the skyrmion core (see Ref. 36 for a detailed calculation). The corrected energy 
is determined by the value p(0) and has the form 



^sk.pin - -gi^Bti T- 



, x 2 / , p (-5/2)' 
(-5/2)/^ „ / *B \ / oz? (-5/2) c dt Q 



R*J \ di 



q T = ±1 . 
(4.10) 



It is assumed here that g 1. 

The skyrmion core radius R* is found by considering the competition between the Zeeman 
energy (4.8) and the energy of Coulomb repulsion. 11 ' 44 Generally speaking, in our case in 
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order to find R* we should include the energy -E s k,pm in the minimization procedure. However, 
this correction only insignificantly influences the result due to the condition R* ^> l B and 
because of the rather strong e-e interaction resulting in the skyrmion formation. Using the 
realistic estimate for the kinematic exchange energy \V\ 2 / A< E c the minimization yields 
the same formula 

E z = £ f(R*/l B ) 2 \n(l 2 B E c /e* z R* 2 ) (4.11) 

as in the case of 'free' skyrmions 44 , where R* 3 = 9n 2 l B e 2 / [64£|«m {E c /e^)\. The number 
of spin-flipped electrons turns out to be rather large 

i / q 2 2\ 2/3 

K= ¥q\M~ b ) ~10-20 (if 4 -0.1K). (4.12) 

We first consider the regime where there is no skyrmions in the clean system but these 
collective excitations could be created due to strong enough kinematic exchange interaction 
between the LL electrons and magnetic impurities. This is the situation where the inequality 
(4.4) is valid, and besides the condition \N — N^\ CiVj is realized, where iVj is the number 
of impurities. The electroneutrality of the system requires that the topological defects are 
created as skyrmion-antiskyrmion pairs. Two impurities are able to create a skyrmion- 
antiskyrmion pair provided the pinning energy E^ pin + E^ pin exceeds the energy increase 
due to the skyrmion-antiskyrmion gap. The latter in accordance with Eq. (4.7) includes the 
Coulomb-exchange part equal to M~ l and twice the Zeeman energy [Eqs. (4.8) and (4.12)]. 
In addition, the energy of a skyrmion and an antiskyrmion pinned by two neighbouring 
magnetic impurities has to be lower than the double energy of a pinned spin exciton. Thus 
the condition Goo < Gi for creation of a pinned skyrmion-antiskyrmion pair can be rewritten 
in the form 

£<£ c ,co, (4.13) 
where the critical value £ Cj00 can be obtained with the help of Eq. (4.10): 

9 + Ft 5/2) (Zcoc) + W2 + E Z A/\V\ 2 = if ^(kco) • (4-14) 

Under the condition (4.13) an impurity acquires the localized magnetic moment K ~ B^ 1 ^ 3 
antiparallel to its own moment and exceeding it (when, e.g. K > 5/2 in the GaAs:Mn case). 

B. Phase diagram of QHF ground state at g^va > 

There are two critical transitions in our problem: first, the global vacuum is destroyed 
when £ becomes less than £ c i and single spin-flip exciton appears (this state may be charac- 
terized as a 'local pinning'); second, the massive pinning of 2DEG spins takes place when £ 
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reaches the value £ Ci00 - However, this scenario is somewhat changed if one takes into account 
finite ratios JVj/iV^. Indeed, up to this point we have supposed that the Zeeman energy e^N^ 
corresponding to the 'global flip' of all 2DEG spins is larger than any contribution to the 
QHF energy due to the magnetic impurities. This global spin-flip actually represents the 
spin configuration treated as the ground state in the previous section. When counted from 
the global vacuum, its energy per one impurity is E s _ { = gifi B B + N^e^/ Ni+\V \ 2 Fq 5 ^ 2 \^)/A. 
Negative E s _{ means that available magnetic impurities are able to polarize completely all 
2DEG electrons even at positive g^oEG- m agreement with the above discussion, one can 
conclude that such a complete polarization takes place when 

#s-f(0 < £min(0 , (4.15) 

where E min = min{0, |VfFi 3/2) (0/A+£|, g^ B B + \V\ 2 F^ l2 \^l A + ^M^ + E^}. 
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FIG. 6: Phase diagram illustrating the reconstruction of the QHF ground state at j 2 * DEG > for 
two cases: isolated impurity Ni/N^ — > (see elucidating legend above the main picture) and finite 
impurity concentration Ni/N^ = 0.01. The calculation was carried out for the Zeeman parameters 
9 = 9il J 'BBA/\V\ 2 = 0.25 and e| = 0.05(g* DEG /0A4)\V\ 2 / A and for the spin exciton dispersion equal 
to £e(q) = 0.05(si* DEG /0.44) + 2£ l-e- q2 / 4 I (q 2 /4) in \V\ 2 / A units. Comments in the figure refer 
to the Ni/N^ = 0.01 case. See text for other details. 



The phase diagram of our system at zero temperature is determined by the interplay 
between Zeeman splitting, Coulomb interaction and kinematic impurity exchange energy, 
and controlled by the impurity concentration. These factors are characterized by the di- 
mensionless parameters g* DEG , £ and Ni/N^. One can construct this diagram by employing 
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the inequalities (4.4), (4.13) and (4.15). The results for both cases of infinitely small and 
finite ratio Ni/N^ are presented in Fig. 6 in the (^* DEG ,() coordinates. We expressed the e-e 
interaction values entering the skyrmion Zeeman energy [Eqs. (4.11) and (4.12)] in terms 
of the parameter M x : E c = M~ l , e 2 ' j 'nl B = M~ 1 {%/-k) 1 / 2 . The phase diagram for the 
Ni/N^ — > case is explicated by the legend above the main picture. The £ ^ci^deg) curve 
in Fig. 6 separates states with unbroken global vacuum (the area above this line) and states 
of 'local pinning' where each impurity is dressed by one bound spin-exciton. The dotted 
line [Cc.oo^deg) curve] separates the state with 'local pinning' and the state of massive spin 
reversal (K ^> 1) determined by the pinned skyrmions (hatched area below this line). 

In the more realistic case of N i /N ( p = 0.01, the curves Cci^deg) and ^coo^deg) formally 
remain the same since the parameter Ni/N^ does not enter Eqs. (4.5) and (4.14). However, 
in this case essential part of the ((^deoO area belongs to states where the 2DEG spins are 
globally polarized in the B direction in spite of positive g* DEG . This area filled by dark-grey 
presents solutions of inequality (4.15). Unbroken global vacuum occupies only the blank 
sector in the upper right corner of the phase diagram. At large £ but fixed |V| 2 /A, the line 
separating the blank and dark-grey sectors tends to g* DEG = 0.088, which corresponds to value 
£z = Ni\V\ 2 /N<pA being the result of the .E s _f(£^oo) =0 equation. At the same time if the 
£^oo limit is realized owing to vanishing V, then both systems of the impurities and of the 
2DEG become independent and at any positive (^dec the global vacuum presents certainly 
the ground state. The light-grey area below the £ ^ci^deg) ^ ne but above the dotted line 
corresponds to the singly spin-flip states with one exciton bound to impurity. The hatched 
light-grey domain below the dotted line corresponds to the state with the localized skyrmions 
created by strong kinematic exchange [Eq. (4.13)]. In our specific case of the Ni/N^ = 0.01 
ratio the dark-grey sector is not contiguous to this skyrmionic region. The total QHF spin 
S z in various states of the phase diagram is indicated in the picture. 

Now we discuss the regime where free skyrmions are already available in the system 
because the number of electrons well deviates from the quantum flux number. Namely, we 
consider that \N — >iVj (although still \N — JV^| <C A^). In this case 'excessive' skyrmions 
may be bound to an impurity. The result depends on the QHF phase. In the globally pinned 
phase (dark-grey area) the binding is impossible since the effective interaction between the 
impurity and the skyrmion is repulsive. In the state of local pinning (light-grey unhatched 
domain) the binding also does not occur. Indeed, the binding energy would be equal to 
E s k,pm (4-10) but due to the condition (4.4) this value is smaller than the spin exciton 
derealization energy -\V\ 2 Fl 3/2) /A. At the global vacuum (blank sector) the binding takes 
place and the binding energy is equal to the pinning energy (4.10). Certainly, the binding 



24 



takes place in the skyrmionic ground state (light-grey hatched sector). However, in contrast 
to the \N— iV^| <Ni case, now all Ni impurities bind skyrmions of the same charge q T , where 
q T = ±1 if correspondingly N^-N^. 

To conclude this section, it is worthy to remind that we have only considered the situation 
where the g* DEG > ground state is realized in the most symmetric phases when the pinned 
spin K is equal to 0, 1 or K ^> oo. As it has been seen, there are only two critical 
parameters £ c i and £ Cj00 in this case. However, one might suppose that transition from 
the local pinning (K — 1) to the skyrmionic phase of massive pinning would proceed more 
smoothly with diminishing parameter £. Namely, below the £ c i(<72deg) curve there should 
be critical value £ = ^(Sjdeg) a ^ which the transition K—l — > K = 2 occurs. This value 
would be the root of equation Gri(£ C 2) = G2{^ C 2)- The next critical point would correspond 
to the K = 2 — > K = 3 transition and so on. This sequence of values £ c i > £ C 2 > ■■■C, c k > ••■ 
where Gk-i(£c^) = Gk^ck) should condense in the vicinity of £ c ,oo- Actually this means 
that the light-grey unhatched domain in Fig. 6 would present not only the singly spin-flip 
2DEG state but a set of states with K = 1,2,3, ... spin-excitons localized at the impurity 
where K is growing with diminishing £. In reality, for a finite e%, K reaches the value given 
by Eq. (4.8) at £ = £ Ci00 - Such a 'stratification' of the light-grey unhatched sector would be 
the only qualitative change of the phase picture of Fig. 6. Quantitative changes would be 
presented by appropriate shifts of the £ Ci00 curve and of the boundary between the dark-grey 
and light-grey areas. However, it is clear that these shifts would not be significant. The 
corresponding crossover parameters £ would at least remain of the same order as the ones 
calculated with the help of Eqs. (4.14) and (4.15). 



V. DISCUSSION 



We have found that the interplay between the kinematic impurity exchange and the 
Coulomb interaction in 2DEG results in the appearance of bound exciton states and in the 
renormalizaton of impurity spin states, including the reconstruction of the QHF ground state 
a t 52DEG > 0. 

Among the available experimental techniques, the inelastic light scattering (ILS) method 
seems to be the most useful method for experimental study of the 2DEG spectra (see Refs. 
7,10,45 and references therein). However, this tool has some special features, and it is 
helpful to discuss our results from this point of view. Let us first consider the (? 2 deg<0 case. 
When measuring the energy from the ground state level Eq , where the impurity has the 
maximum spin projection [see Eq. (3.10)], one obtains ten levels of the localized excitations 



25 




FIG. 7: The case of <7 2D eg < 0. Energies of the excitations measured from the ground state 
are plotted in units of |V| 2 /A as functions of £ = A/M X \V\ 2 . Filled areas show energies of the 
delocalized spin excitons. See text for details. 

AE related to the spin changes SS Z = 0, 1, ...5 

AE 5Sz ,t = 9i» B B5S z + (|Vf /A) [i?W(fl _ F Q ( - 5/2) (0] (t = 0,x), (5.1) 

where s = 5S Z — 5/2, and the index t labels the type of the excited state [in Eq. (5.1) it is 
considered that F^ 2 ^ = 0]. Within the scope of the experiment where only the \5S Z \ < 2 
excitations seem to be observable as ILS peaks, we plot in Fig. 7. these five levels calculated 
with the help of Eq. (3.9) for the parameters g = 0.25 and ez = 0.05|V| 2 /A. This calculation 



is done for the sake of demonstration with the function S q = 2M X 



l-e-^ 4 / (g 2 /4) 



and 



the fitting parameter M x used to describe the spin-wave dispersion. In the available wide 
quantum wells the inverse spin-exciton mass is relatively small. 46 Hence the values £ < 1 
seem to be experimentally relevant, and the evolution of non-equidistant excitations AEss z , x 
as a function of £ (and therefore of B) should be observable in this interval. 

The non-localized states discussed in Section III C actually present a transformation of 
the x-type excitations when the spin exciton is detached from the impurity and falls in the 
spin-wave continuum. The bottoms of continuous bands are shown as filled areas in Fig. 7. 

(<5 S z —— ) 

The band edges are higher than the AEss ZlX curves by the quantity —(\V\ 2 /A)F K z 2 (£) 
[see Eqs. (3.8) and (3.10)], and therefore the latter may be treated as spin-exciton binding 
energy. However, it seems to be difficult to observe these states in the ILS spectra because 
of comparatively small oscillator strengths, specifically, due to divergence of the envelope 
function f(q). 

Similar ILS picture should also take place for g* DEG > in the phase of the 2DEG global 
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pinning (dark-grey area in Fig. 6). In the skyrmionic phase (light-grey hatched sector) 
there are intra-impurity ILS transitions determined by Eq. (5.1). Besides, two new types 
of resonances are expected: the first is the skyrmion derealization with 5S Z = and with 
excitation energy equal to E skpin (4.10); another one is transition 5S Z = — 1 where the 
delocalized skyrmion leaves the impurity with the bound spin exciton. In the latter case the 
transition energy is E skpin — \V\ 2 F^^ +e^- 

In the global vacuum and local pinning states (blank and light-grey unhatched domains) 
the ILS spectrum is determined by transitions between levels (4.1) and (4.2), so that, e.g., the 
ILS transitions to the localized states from the global vacuum are determined by the energies 
AE SSz = e* z + giH B B(5S z — l) + IV^F^ 2 ^ 3 ^ /A an d correspond to nonzero spin change 
5S Z = 1,2, (3, 4, 5). At the same quantum numbers SS Z there should also be resonance features 
related to the impurity spin rotation, which cost the energy AEss z ,res — 9iHbB8S z . These 
resonances are in fact transitions to the continuous spectrum which should be noticeable on 
the background of free spin waves contribution (e z < AEkE^) [see the comment above Eq. 
(4. 2)]. 47 The ILS spectrum of excitations from the local-pinning ground state is presented: 
first, by the 8S Z = — 1 transition to the global vacuum (this energy is equal to AE_i = 
— |y| 2 Fx 3 ^/A — second, by the 5S Z = 1,2, (3,4) transitions to the localized spin-flip 
states with energies AE 5Sz = gifi B B5S z + \V\ 2 Fi 3/2_55j) -F^ ,2) /A; and third, by the 
5S Z = 0,1,2,(3,4) transitions to the resonance states in the continuous spectrum with 
transition energies AE SSztTCS = gifi B B(5S z + l) - |V| 2 Fi 3/2) /A - £ * z . 

Finally, we note that currently there are several possibilities for the experimental study of 
skyrmion-like textures (e.g., see Refs. 13,14,17 and 48). However, for any method one of the 
most serious obstacles impeding observation of spin-flip phases is the very narrow interval in 
the vicinity of the s^deg — factor where the spin-flip reconstruction of the ground state or 
skyrmion-like excitations are possible. From this point of view the minor magnetic doping 
would become an additional fine tuning tool allowing to change the balance between E z and 
Eq and to influence the skyrmion formation. 
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APPENDIX A: EXCITONIC REPRESENTATION 

The excitonic representation is a convenient tool for a description of electron-hole states 
in a 2DEG multiply degenerate in m. When acting on the vacuum state |vac) (in our case 
this vacuum is defined in Subsection II A), the exciton-creation operators form a system 
of basis states diagonalizing the Hamiltonian including a considerable part of the Coulomb 
interaction. Due to translational invariance of a clean 2DEG these exciton states are classified 
by the 2D momentum q and the degeneracy turns out to be lifted. The exciton operators 
for a single LL were first introduced in Refs. 6. The commutation rules for the same case of 
single LL were found in Ref. 49 (see also Ref. 9 and references therein). 

Unlike previous papers, where the ER technique was developed for the Landau gauge, 
we use the symmetric gauge for bare one-electron states. In the symmetric gauge the spin- 
exciton creation operator is expressed in terms of the a m and b m Fermi operators, 

N+-1 

Cq = E h* mm ,(q)bi_ n a m ^ , (Al) 

m,m'=0 

(cf. the definition based on the Landau gauge 6 > 9 > 18 > 36 ). I n this expression 

hmk(q) = (m\/k\) 1/2 (q^ k - m L k m - m (q 2 /2)e-i 2 / 4 (A2) 

are the building block functions used in the ER technique, g_ = %qe~ %lp j y/2 = i(q x —iq y ) / \/2, 
and L^ m are the Laguerre polynomials. Here and below all lengths are measured in the 
magnetic length Ib = 1 units. The spin-exciton states are orthogonal and normalized, 

(vac|Q qi Q q2 |vac) = <J qi>qa . (A3) 

The operators (Al) together with the intra-sublevel operators 

A * = j\T E hmm'^Wm-nam'-n, and ^ = Yl k ™™' i^m-rfim' -n ■ (A4) 

^ m,m'=0 ^ m,m'=0 

form a closed Lie algebra. In order to check it we first obtain the commutation relations 
[Ql 4J = N~ 1/2 h* m+n , k (<l)bt n , [Ql M = -N; 1/2 K, m+ n(^W-n, (A5) 

fc=0 k=0 
N*-l . Nj,-1 



[A,, am] = -— >W,fc(<j) a fc-ii, [B q , b m ] = — /i fc , m +„(q)&I_ n , (A6) 

and 

[Qq, 4J = [Qq, 6m] = [A,, &m] = [A,, *4] = [#q, «m] = [Bq, 4j = . (A7) 
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As a result we see that operators (Al) and (A4) really form a closed algebra with the 
commutation relations 9 ' 18 ' 49 

\0 0+ 1 — P *(qi xc l2) z /2 A _ „-i(qiXq 2 )z/2» 

e- i(qiXq2) ^ /2 [<, Q q2 ] = -e^*<^/ 2 [£ qi , Q q j = -A^Q^, 
«^ q J = -^sin[( qi xq 2 ) z /2]4 ( A8 ) 



qi'^q 2 J _ jV^ H2^/-"J^qi-Ki2) 

[^ qi ^qJ = ~ sin [( qi x q 2 ) 2 /2]3 qi ^ 2 . 
Acting on the vacuum state the intra-sublevel operators result in 

-4 q |vac) = 5 qj o, and i3 q |vac)=0. (A9) 
The excitonic basis Q q |vac) determine the set of eigenstates of a clean 2DEG, 

H[ s) + if M ), Q q ] |vac) = (e z + 5,) Q q |vac) . (A10) 
Here £ g stands for the Coulomb energy of the free spin wave defined by the equation 3-5 

£ « = wH W ^ i l - e * (PXq)z ] = / dPwM ^{P)} 2 [1 " JoiPQ)] , (AH) 

Jo(pq) is the Bessel function. 

The Coulomb vertices in the Hamiltonian (2.12), (2.13) are given by the equations 

W ss {q) = v ss (q)[h nn (q)} 2 , W sd (q) = v sd (q)h nn (q), (A12) 

where 27rt> ss (q) and 2irv sd (q) are the 2D Fourier transforms of the average s-s and s-d 
interaction potentials. One can present them as 1 

v ss (q) = ^- q jj dz 1 dz 2 e-^-^\az 1 )\ 2 \az 2 )\ 2 , 

v sd {q) = ^~ q J dze-^\\C{z)\ 2 . (A13) 

(The impurity cite is assumed to be at the point = {0, 0, z d } .) 

APPENDIX B: SPIN OPERATORS 



Bound spin excitons are characterized by the spin numbers S z and S . The corresponding 
operators have the form 

S z = Si s) + Si d) , (Bl) 

where 

= ^ (A - Bo) ; (S«) 2 = N^Qo + (S^) 2 + (B2) 
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and 



j(n T -ni); (S (d) ) 



23 3 

n T + - -n T nj. 



(B3) 



2 v-, > v' / 4 V "' ' " l/ 2 

are the spin operators for 2DEG (in the excitonic representation) and for magnetic impurity 
(in terms of the single-orbital model), respectively. The total squared spin operator for the 
system is defined as 

2 ,,„/,, , \ , ,\ 2 



S 2 = (s«) + + < 2 (Qt c t q + c { CtQo ) + (s^y 

The operator S z commutes with the Hamiltonian (2.9), while for S 2 one has 



(B4) 



S 2 ,H\ =iV] /2 (^/i BJ B-£ Z )c}c T Q 



+ 



(ft - ft)V [c{ (n T 6 + iV; /2 Q ao) - c{ (^a + ^^Q^^ 



(B5) 



H.c. 



The difference between the g factors of the magnetic impurity and the host QHF and the 
difference between the projection factors ft and ft measure the spin non-conservation. 
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